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ABSTRACT 

A reexamination of the effects of non-zero degeneracies on Big Bang Nucleosynthesis is 
made. As previously noted, non-trivial alterations of the standard model conclusions can 
be induced only if excess lepton numbers comparable to photon number densities T) r , 
are assumed (where t) r ~ 3 x 10 9 t 7 6 ). Furthermore, the required lepton number densities 
(D, - 77r) must be different for v t than for v ^ and u T . It is shown that this loophole in the 
standard model of nucleosynthesis is robust and will not vanish as abundance and reaction 
rate determinations improve. However, it is also argued that theoretically | L t | ~ | \ 

~| L t | ~ r/i, <C rj T which would preclude this loophole in standard unified models. 
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The baxyon number of the Universe is a quantity of fundamental interest. Baxyon 
number violation is a process of considerable importance in both cosmology and parti- 
cle physics. In many grand unification models, as well as in recent studies of the elec- 
troweak model at high temperatures, violation of B is accompanied by violation of lepton 
number L, generally with conservation of some linear combination of the two ( B — L 
for instance, in minimal SU5). In cosmology, the baryon to photon ratio, though small 
(77 = (ng — ng)/n T = 0(1O -10 )), has a major influence on the primordial abundances of 
the light elements; may determine when the Universe became matter-dominated; and is 
an important parameter in theories of galaxy formation and dark matter determinations. 
Clearly it is a quantity whose value we would like to know. Fortunately, standard big bang 
nucleosynthesis can provide us with relatively tight constraints, 7710 = ^7/ 10 ~ 10 = 2. 8-4.0 
[1]. If lepton number violation is of the same order as baryon number violation then it has 
a negligible effect on nucleosynthesis. However, the question remains, can nucleosynthesis 
give any constraints on the violation of lepton number? To answer this requires a rela- 
tively straightforward adjustment to the usual calculation of the primordial abundances of 
the light elements. Although large lepton asymmetries can cause difficulties with models 
for galaxy formation [2], they do not preclude the possibility of galaxy formation when 
additional physics (e.g. late time phase transitions [3]) is taken into account. 

Since we know from redshift-luminosity measurements that the total density of the 
Universe is not much more than critical we can infer a limit on baryon number from 

n B hl0~ 3 = 3.53 x 10 7 77 (1) 

where the Hubble parameter Ho = holOO km s _1 Mpc -1 and 6 is the microwave back- 
ground temperature in units of 2.75 K. If we assume that ft = 1 (fttotal> n °t ftfl) then age 
of the Universe arguments constrain ft/ig^ 0.25 [4] and thus 7710 <; 70. More generally, we 
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can safely assume that ^ 1 which implies 7710^ 280. In addition, the upper limit on 

the fractional excess of charge in the universe (n c — is of order 10 18 [5] so we can 

assume that the lepton number due to electrons and positrons is small. Any additional 
asymmetry in the electron family must therefore be in the form of neutrino degeneracy. 
The approach we take then is to allow the three known neutrinos to have chemical poten- 
tials (& = iXyJTv) up to the value 54, at which point the energy density of that species 
reaches the critical value for closure, p c = ZHq/8kG (see equation 3 below). 

The exploration of neutrino chemical potentials for big bang nucleosynthesis has been 
performed on several occasions in the last few years [6-10], however a number of the 
relevant reaction rates in the calculation have been updated recently as a result of improved 
measurements. In particular, since the recent calculation of Terasawa and Sato [8] an 
improved measurement of the neutron half-life [11] dramatically narrows the uncertainty in 
the n<->p rates [1,12]. In addition the rates D(d,n) 3 He, D(d,p)T, 4 He(t,7) 7 Li, 7 Li(p,a) 4 He 
and 7 Be(n,p) 7 Li have changed [13]. The purpose of this paper is to investigate the effect 
of these improved input parameters on the baryon and lepton asymmetries of the Universe 
to accurately assess the current situation. 

In addition to improvements in the numerical computation of primordial abundances, 
the comparison between the predicted abundances and the observational determinations 
has become more stringent. For example the older bounds on the 4 He abundance Y p < 
0.25 or 0.254 have now been improved to Y p = 0.23 ± 0.01, and the bound on the lithium 
abundance which has ranged from 7 Li/H < (3 — 10) x 10“ 10 is now down to 7 Li/H<; 1.4 x 
10” 10 [1]. In what follows we will examine the effect of these changes on the neutrino 
degeneracy loophole. We will conclude that the situation is not qualitatively different from 
the previous conclusions, despite more stringent input data. Furthermore, we will show 
that loopholes in the standard big bang nucleosynthesis conclusions about baryon density 
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will always persist if one allows the introduction of different but specific combinations of 
values for and and/or £ r . While the values required will seem unphysically large, as 
we will discuss, this loophole will nonetheless remain. 

Our calculation is based on the code of Wagoner et al [14], with the n<-+p rates 
calculated as described in [1]. We also use the most recent reaction rates of Caughlan and 
Fowler [13] and the recent data on the neutron mean-life. The observational limits that 
we use on the light elements are those found in ref. [1]: 

0.22 < Y p < 0.24 (2a) 



where Y p is the mass fraction of 4 He and X refers to number density relative to hydrogen. 
The most recent calculations [1] combined with the above observational bounds indicate 
that the simplest version, ie. the standard model of nucleosynthesis, is consistent with 
these observations. The consistency occurs when 2.8 < t]\q < 3.3(4.0). (The higher value 
up to r ] io = 4.0 is allowed when the uncertainties in key 7 Li rates are taken into account, 
see eg. ref. [1,15].) Through the years, there have been numerous attempts at altering the 
conclusions of the standard model, by (usually) complicating the input. For example early 
claims of allowing fig = 1 because of inhomogeneities or decaying particles have largely 
been dispelled. The nagging alternative of including neutrino degeneracy is the one we 
discuss here. 

Introducing neutrino degeneracy has two effects: (1) weak reaction rates are altered by 
the change in the electron neutrino distribution function, thus changing the equilibrium 
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ratio of neutron, to proton densities, and (2) the energy density of neutrinos increases, 
speeding up the expansion of the universe. Thus £ e has two effects, whereas and £ r only 
affect the expansion rate. Furthermore, this means that and £ r axe interchangeable as 
far as their effects are concerned. 

The energy density of neutrinos (and antineutrinos) can be evaluated analytically: 


«» = £ T ^) 4 






£« = Mi *ilTu 


(3) 


where we assume the tau neutrino is effectively massless for T ~ 1 MeV. For a single 
neutrino species, closure density is achieved for £ = 53.8 (for a Hubble parameter of 
100 km s -1 Mpc -1 ). The lepton number of the universe is then 


Li 


n(vj) - n(vj) 

»( 7) 


1 

12C(3) 



(*■ 2 6+e, 3 ) 


(4) 


(where £ is the Riemann zeta function and £(3) ~ 1.202). If lepton number is conserved 
then & is a constant, except during e + -e - annihilation. 

Since the chemical potentials for fi and r neutrinos appear only in the expression 
for the energy density, introducing £ r is clearly equivalent to introducing additional 
neutrino flavors. We can write 


w^r = 3 + ^^U?+{?/2^). (5) 

where N ue ff parameterizes the neutrinos in terms of massless neutrino species. Following 
the Z° results from LEP we will naturally assume that there are three neutrino species 
and that the tau neutrino is light. Note that although ~ 50 corresponds to N ue ff ~ 10 5 , 
one could equivalently add 10 5 fermion degrees of freedom which do not couple to the Z° 
and hence avoid the LEP limit. 
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The effect of neutrino degeneracy on the helium abundance can be readily understood: 
Increasing or £ r (with either sign) raises the expansion rate before nucleosynthesis. This 
leads to a higher freeze-out temperature for weak interactions and hence an increased yield 
of 4 He. (This of course, is the same argument as for an increased number of neutrino 
flavors.) The dramatic increase in the expansion rate (due to the equivalent of ~ 10 5 
neutrinos) is allowable when one notes that although raising also affects the expansion 
rate (increasing it), in addition it changes the weak reaction rates. The neutron to proton 
density ratio (in equilibrium) is given by 

n/p = exp(-A m/T - £«) (6) 

where Am = m n — m p — 1.29 MeV, so increasing £ e leads to a smaller value of n/p when 
the weak rates freeze out and hence a smaller yield of 4 He. Hence, for fixed 77 the two 
quantities £ e and can be played off against each other without grossly affecting the 

4 He abundance. Changing the sign of has the opposite effect on the weak reaction 
rates, but gives the same contribution to the energy density. As 77 increases so does 4 He 
production as the nuclear reactions producing 4 He become more efficient relative to the 
expansion rate. Increasing the expansion rate (by increasing £ c ,*i,r) and decreasing ( n/p)? 
(by increasing £ e ) can readily compensate for an increase in 77. 

D+ 3 He and 7 Li production are far less sensitive to n/p at freeze-out as their abundance 
is primarily determined by competition between nuclear reaction rates and the expansion 
rate. The longer the nuclear rates are in equilibrium, the more D and 3 He are destroyed. 
For 7710^ 3 (<; 3) the production of mass 7 nuclei ( 7 Be, 7 Li) increases (decreases) with 
increasing 77. Thus an increase in £ e ,/i,r can always compensate for an increase in 77. At 
fixed 77 relatively large increases in £ e (driving n/p 0 exponentially) are necessary to 
compensate increases in In the case of D+ 3 He, the increase in £ e shuts down D 
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and 3 He production. For tjio^ 3 the increase in shuts off 7 Be(n,p), the destruction 
channel for mass 7, and results in more 7 Li. For f/io 3 an increase in £ e results in less 
7 Li production (via 4 He(t, 7 )) and greater destruction (via 7 Li(p,a)) — the net result being 
a decrease in 7 Li. 

In figures la-d we show the regions of the plane (taking £ r = £ M ) allowed by 

the observational limits for rj io = 2.8, 3.3, 10, 280 and for r n = 889.6 sec. The curves in 
the figures are iso-abundance curves at the observational bounds given in (2a— 2d). The 
allowed region is shown in bold. For rji o < 2.7 the region in which the 4 He limits are 
satisfied has X j + -X3}j e > 10 -4 . 

At 7710 = 2.7 we axe able to meet the constraints for ^ = 0 and £ e = 0.1. As 

we increase 77 the yields of deuterium and 3 He begin to drop and at 7710 = 2.8 we are able 
to have £ e = = 0 (fig. la). (If we increase sufficiently (to ~ 25) the increased 

expansion rate reduces the time available for nucleosynthesis and brings Y p back down to 
~ 0.24 [16]. However the D and 3 He abundances continue to rise, so this region is still 
ruled out. This corresponds to ~ 18000 neutrino flavors [17].) 7710 = 2.8 corresponds to the 
lower limit of ref. [1]. Here the limits on the chemical potentials are |^| < 0.6 (if = £ r ) 
and —0.02 < £ e < 0.1. The lepton numbers are then constrained by \L p \ < 0.15 and 
—0.0050 < L e < 0.025. As 77 rises further the deuterium and 3 He abundances continue to 
drop, while Y p begins to increase. (Raising 77 causes nucleosynthesis to begin earlier, giving 
a larger value of n/p at the onset, and hence increasing Yp.) The limits on + -^ 7 Be a l so 
begin to move to higher values of the chemical potentials, leaving the allowed region around 
0£ £ e & 0.1, 0<; £p<; 2 for 7710 = 3.3 (fig. lb). This is the highest value of 77 for which the 
origin falls within the allowed region. As we continue to raise 77 the observationally allowed 
region moves to higher values of £ e and £ M . When we reach 7710 = 280 (where Clg = 1 
for Hq = 100 km s - 1 Mpc -1 ) we find we require £ e « 1.6, ~ 40 in order to satisfy the 
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observational constraints. 


In fig. 2 we show the same contours as a function of 7710 and f e for ^ = f r = 0. Here 
it is clear that for f e = 0 the limits on the baryon to photon ratio axe 2.8 < 7710 < 3.3 in 
agreement with ref. [1]. However a value of f e between —0.01 and +0.09 is still permissible. 

We have seen therefore that the standard nucleosynthesis bound 7710 < 3. 3(4.0) can be 
bypassed by introducing two new parameters and in fact allow for fls = 1. It is amusing 
to note that in addition to the limit on 7710 , the parameters f e and are sufficient 

for obliterating the cosmological bound on the number of neutrino flavors. Indeed by 
choosing = f r = 0 and f e at its limiting value of ~ 1.6, the cosmological bound becomes 
N v £ 1.8 x 10 5 . A further twist arises if for some reason one were compelled to take 
£/*,r ~ 0(f e ). In that event, an extra source of energy density comparable to that given 
by the bound on jV„, N„ ~ O(10 5 ), is needed to account for the necessary expansion. 

Despite the baroque nature of the preceeding discussion, the nucleosynthesis loophole 
based on non-vanishing neutrino chemical potentials does not lead to any direct inconsis- 
tencies, per se. Furthermore it is one that we do not forsee as disappearing, either because 
of improved cross-section measurements or abundance determinations. This can be seen by 
an examination of figs. la-d. The closest observational bound that could close the loophole 
is the upper bound on (D+ 3 He)/H. For (D+ 3 He)/H^ 8 the window in the 7710 = 280 plot 
closes, but then the standard model is also ruled out for any value of 7710 . 

We do not however wish to leave the reader with the impression that we axe advocat- 
ing this route for achieving a large value for Q b , or that one should ignore the cosmological 
bound on N v . In principle it may be possible to have Li (if,) be non-zero, but any realistic 
theory to date would only predict Li ~ 0( 77 ). It would also be peculiar to have L ^ or L T 
much different from L e as would be required to have nucleosynthesis go through the loop- 
hole [18]. Furthermore, unless the baryon asymmetry is generated very late, electroweak 
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baryon number violation which conserves B — L would require that the total lepton asym- 
metry L = L^+Lp+Lr =r}(L = riL/n 1 = ns/n 7 ) [19]. Hence, degenerate nucleosynthesis 
which requires L e ,X M)T ^ tj and | i r [ > L e would necessitate a very special cancellation 
to achieve L = 77. Thus while we cannot categorically eliminate large neutrino degeneracies 
as a loophole for cosmological nucleosynthesis it is clear that they are exceedingly unlikely 
on general theoretical grounds. 
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FIGURE CAPTIONS 

la. Limiting contours in the plane for 7710 = 2.8, £ r = £*» and r n = 889.6 sec. The 

limits shown are those of equations 2. The region of the plane allowed by all limits is 
shown in bold. 

lb. As fig. la but with 7710 = 3.3. 

lc. As fig. la but with 7710 = 10. 

l d. As fig. la but with 7710 = 280 and on a larger scale. 

2. The same limits in the £e -7 7io plane with = 0. 
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